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Abstract. In recent years we've seen the birth of a ncw field known as hamiltonian complezity lying at the 
crossroads between computer science and theoretical physics. Hamiltonian complexity is directly concerned 
with the question: how hard is it to simulate a physical system? Hero I rcvicw the foundational results, 
guiding problems, and futurc dircctions of this emergent field. 



1. Introduction 

In the natural sciences coniplex behaviour is usually exhibited by large numbers of interacting particles 
nioving according to simple underlying physical laws. Many body theory is the area of physics aimed at 
explaining these systems, a task oftcn rcquiring the study of equations comprising large numbers of variables. 
Problems involving many variables like these are a recurring theme throughout the natural sciences and 
there are several approaches to undcrstanding thcir propcrties. The program that concerns us here is that of 
computational complexity theory, which is the branch of computer science and mathematics focussed on the 
task of classifying thcir computational difficulty. 

There are striking parallcls between the ccntral objccts of study in many body theory and computational 
complexity theory. For cxample, the hamiltonians of many body physics bcar some superficial resemblance 
to the propositional formulae studied in complexity theory. It was quickly realised that this analogy runs 
deeper: the equilibrium states of many physical systems directly correspond to satisfying assignments of 
variables in certain families of propositional formulae. FoUowing this it made sense to apply the apparatus of 
computational complexity theory, aimed at answering the qucstion "how hard is it to decide a propositional 
formula?" , to the study of many body physical systems. Although these analogies were explored in some detail 
during the 80s and 90s it was only with the advent of quantum Information theory that we have witnessed 
an intensified cross-fcrtilisation between computational complcxity theory and many body physics, leading to 
the birth of hamiltonian complexity. 

The ccntral objectivc of hamiltonian complexity is thus to answer the question: "/low hard is it to simulate 
a physical systeml" Computational complcxity theory provides us with the apparatus to quantify what is 
meant by "how hard" and many body theory, in particular, condensed matter physics, provides us with a 
rich family of strongly corrclated systems to study. 

Hamiltonian complexity is interesting prccisely because there are families of systems that could be, in 
principle, engineered and whose physical properties are hard to simulate. The canonical example here is due 
to Barahona |15j who described an Ising model in a nonuniform magnetic field whose ground-state energy 
is hard to calculate. By "hard" we mean that the computational resources required to estimate it grow too 
quickly with the system size, so that no computer can carry out the task in a reasonable time. Further, 
the origin of this hardness isn't just a lack of imagination: no matter how clever you are there is no way 
to efiiciently extract predictions about these system|3- This peculiar situation is initially disturbing: if one 
cannot solve the model in a reasonable time, even with access to a supercomputcr, then how can the model 
be said to have predictive content? Could it be that there exist systems whose behaviour is simply too hard 
too predict? This is reminiscent of chaos theory, but the origin of the lack of predictability is rather different. 

So far this discussion applies only to classical systems. However, it was realised that guantum computers 
would be capable of simulation tasks excccding thcir classical countcrparts. Indeed, a quantum computer can 
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'^This is because merely the existence of an efficient method to compute the ground-state energy of Barahona's model would 
serve as a counterexample to the P 7^ NP conjecture. I'U say more about this later. Obviously, for this review, I am adopting 
the position that this conjecture is truc. If it were falsc then this would have some spectacular conscquonccs for the natural 
sciences... 



1 



2 



TOBIAS J. OSBORNE 



simulate the real-time dynamics of naturally occurring many body systems, both bosonic and fermionic |76j . 
With this obscrvation came the hope that quantuni computers would also be ablc to simulate the eguilibrium 
dynamics of, e. g., Barahona's model, and much effort was expended on the task of developing such a quantum 
simulation algorithm. Unfortunately - or fortunately, depending on your point of view - it became increasingly 
evident that such quantum algorithms couldn't always suceeed [T]. Concomitant with these investigations 
has been a parallcl effort at classifying the guantum complexity of computational problems, eulminating in 
the construction of quantum models which at first sight appear to be simulatable with a quantum computer, 
yet whose properties are, in a specifie sense FH deseribe later, "even more" intractable than their elassical 
cousins. 

The puzzling obscrvation that there are systems whose equilibrium properties are hard to simulate has 
a relatively simple resolution: it turns out that the assumption that the system ends up in an equilibrium 
state modelled by the Gibbs distribution is simply wrong. This isn't really surprising, and is familiar in 
the context of spin glasses |16j : when a system has access to a large number of metastable configurations 
it is overwhelmingly to end up in a metastable state rather than the unique lowest-energy state. Indeed, 
Barahona's model, and its cousins, exhibit precisely this property. Thus, if you were to engineer a "difficult" 
system such as Barahona's, it would take an extremely long time for it to equilibrate to its ground-state, and 
any elassical or quantum simulation algorithm for its ground-state properties would get stuek. Thus we've 
learnt that the manifold of states accessible by the natural dynamics of a physical system is a tiny subset of 
the full configuration space of the system. By understanding this corner of configuration space better we can 
design better variational classes of wavefunctions to best parametrise naturally occurring states. 

The outline of this review is as foUows. We begin in ij2]with a brief diseussion of computational complexity. 
In ij3]wc discuss the ccntral question of hamiltonian complexity, the simulation problem, and review its origins 
and motivations. In i]4]we review a large body of literature focussed on the presentation of physical systems 
that are hard to simulate. Finally, in ij5l we review the literature focussed on the converse situation of 
understanding what physical systems can be efhciently simulated. 

2. COMPLEXITY THEORY 

In computational complexity (see, e. g., [81] or [10] for a textbook treatment) one categorises the difficulty 
of computational problems by organising them into complexity classes. To understand which class a problem 
belongs to you need to study how many basic computational steps (e. g., addition, subtraction) are required 
to solve the problem for larger and larger inputs. It is typical that the larger the input to the problem is 
the more computational steps are required to solve it. For example, consider the problem of sorting a list 
of n numbers in descending order. As the input size n gets larger it takes, in the worst case, no more than 
0{nlog{n)) steps to do this. Anothcr example: consider the problem of multiplying two n x n matrices. 
This takes, in the worst case, no more than 0(n^+'^), e < 1, operations to complete [29]. In computational 
complexity one doesn't care so much about the exponent on the n, rather, that there is some finite exponent 
a < cxD so that, in the worst case, no more than 0(n") operations are required to solve the problem. 

The class of all decision problems - i. e., problems with a "yes" or "no" answer - solvable using a number 
of basic computational steps scaling no worse than a polynomial in the input size n is called P. The exact 
polynomial rate needed to solve a problem in P may differ from computer architecture to architecture (the 
actual definition of P doesn't depend exactly on the definition of computer), but the fact that the rate is 
polynomial doesn't. This is an admittedly coarse classification because for all practical purposes a problem 
requiring n^^^'^ steps to solve is intractable for large instances. Interestingly, for many problems in P the 
exponent seems to be low. 

Another fundamental complexity class is NP which is, roughly speaking, the class of decision problems 
whose solution may be checked using a polynomial number of basic computational steps. At first sight 
this innocuous definition doesn't hint at any hidden intractability. However, it turns out that N P contains 
extremely complex problems. Firstly, it is worth noting that P C NP since the solution to any problem in P 
can be checked using polynomial steps (just solve it yourself). To convince yourself that N P might contain 
more problems than P consider the foUowing situation. Imagine you have some funetion f oi n bit values 
xi, X2, . . . , Xn, Xj G {0, 1}, and you want to find its minimum value. (One can easily convert a problem such 
as minimum finding into a decision problem by solving the related problem: is the minimum of / less that 
c? If you can solve this problem then you can use a binary search to find the minimum value up to some 
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precision.) It is convenient to imagine that the bit valucs x = (2:1, X2, ■ ■ ■ ■, Xn) represent thc orientations of a 
collection of n spins and that / represents an energy function of x. Now, as long as it is possiblc to compute 
/(x) using only a polynomial numbcr of operations, onc can see that this problem of finding thc minimum 
is in N P. (That is, if someone gave you the configuration which is claimed to achieve a value of / less than 
some c then it is easy to check that this is the case.) But is it in P? 

There is one obvious strategy to find a value of / less than some desired c: simply search over the space 
of all bit strings x. Howevcr, this task is prohibitively expcnsivc (it takes 0(2") stcps in thc input size n). 
Perhaps there is a better strategy? Unfortunatcly not, unlcss / has some special structurc. For cxample, if 
/ is compriscd of a scqucnce of random intcractions - cach coupling collcctions of many spins - the system 
will be frustrated and / will posscss many local minima, ruining any local search heuristic. 

The canonical example of such a frustrated / is provided by the 3SAT problem: here / is given by a 
Boolean expression writtcn as a conjunction (i. e., the "OR", denoted V) of a collection of clauses Cj, each 
of which is the disjunction (i. e., thc "AND", denoted A) of thrce literals (i. e., a Boolean variable Xj or its 
negation Xj). An example of an allowed clause is, e. g, Cj{x) = Xj-^ V a;^^ V Xj^. There are 8 possible choices 
for which variables can be subjected to flips and all are considered to be possible. You should think of which 
variables are actually subjected to flips in each clause, and which variables appear in each clause, as being 
choscn by an adversary. The function / corresponds to a (classical) hamiltonian for a system of n spins with 
general 3-spin intcractions via 

m m 

(2.1) /= /\C,(x)^F(x) = -^/i,(x), 

where each of literals Xj are interpreted as numbers via TRUE = 1 and FALSE = O and the intcraction 
terms h j are found by replacing V with the Standard product operation, e. g., /ij(x) = Xj^Xj^Xj^ i — !■ Cj(x) = 
Xj^ VXj-2 Vxj3 . A configuration of zeros and ones achieving the minimum value of H corresponds to a satisfying 
assignment of / where all of the clauses evaluate to TRUE. If zero is regarded as a spin up and one is regarded 
as a spin down then H can be thought of as the hamiltonian for a (possibly frustrated) spin system. 

It is not impossible that there is a general-purpose efScient algorithm to find the minimum value of H, but 
this would have fabulous consequences because of a fundamental result in computational complexity theory 
known as the Cook-Levin theorem: the 3SAT problem is complete for the class N P. The theorem says that if 
you had an algorithm which could solve 3SAT using a polynomial number of steps, in the worst case, then you 
could solve any other problem in NP efficiently. Any such problem in NP with the property that an efRcient 
algorithm for it would yield an efficient algorithm for all of NP is known as an NP- complete problem. A 
problem is known as NP-hard if it as at least as hard as the hardest problcms in NP. There is no requirement 
that such problems be in N P; indeed, they may not even be decision problems. The class NP is now known 
to contain a vast number of problems, many of which have considerable practical relevance. However, since 
no efficient solution algorithm has ever been found, nor are there scrious proposals, the possibility of P = N P 
is considered to be extraordinarily unlikely. 

Interestingly, there are some surprising examples of systems that appear, at first sight, to be intractable, yet 
which admit efficient algorithms for their physical properties. A good case in point is that of a ferromagnetic 
Ising model - where the spin-spin intcractions are always positive - with arbitrary interaction graph. Thanks 
to Barahona's example (which pertains to an Ising model with both antiferromagnetic and ferromagnetic 
intcractions), one might be tempted to assume that the equilibrium properties of this model cannot be 
simulated efficiently. However, there is a clever polynomial-timc algorithm which samples configurations 
according to its Gibbs distribution [51], to a required precision. 

With the advent of quantum compution the landscape of computational complcxity theory has been fun- 
damentally extended: it is natural to introduce a new class of problems extcnding the classical clasfH P, which 
are easy for a quantum computer to decide. This natural class is known as BQP and contains several problems 
thought to lie outside of P, e. g., factoring. The classification of the computational complexity of problems in 
the presence of quantum computers is known as quantum complezity theory. 

There is also a natural class of problems known as QMA, analogous to N P. To understand this generalisa- 
tion it is convenient to view N P as a game played between two players: (M) a malevolent Merlin, with infinitc 

"^Strictly speaking, the class BQP extends the class BPP, which are the decision problems that can be decidcd with a boundcd 
probability using a polynomial number of arithmetic operations and coin flips. 
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computational power; and (A) Arthur, limited to using a Standard classical computcr with polynomial com- 
putational power. In this game Arthur asks Merlin questions, such as "does / havc a satisfying assignment?" 
and Merlin accompanies the answcr with a proof x. which Arthur can usc to vcrify thc answcr. In the case 
where / doesn't have a satisfying assignment there is no way Merhn can ever fool Arthur that it does. The 
class NP is then the class of ah decision problems which can be decided via this game. 

The quantum generahsation of NP is defined along the same hnes, except that Arthur is ahcwed access 
to a guantum computer with polynomial computational power and Merlin is allowcd to scnd quantum states 
as a proof. Arthur asks the same questions as before, but now Merlin is allowcd to scnd a quantum state to 
convince Arthur that a given problem has a "yes" answcr: Arthur accepts that a given problem has a "yes" 
answcr if, after running a quantum verification algorithm on Merlin's quantum state, he measurcs a 1 and 
rejects if he measurcs a 0. We say a decision problem is in QMA if Arthur always accepts a "yes" instance 
with probability greater than 2/3 and always rejects a "no" instance with probability less than 1/3. 

A remarkable result of Kitaev (discusscd bclow) cstablishes a Cook-Lcvin-likc theorem for QMA, and 
provides a complete problem for it. Thc main intuition to take away here is that the class QMA contains 
"harder" problems than N P. 

3. The simulation problem 

Typically the way physics has done is as follows: (1) formulate a theoretical model; (2) (approximately) 
solve the model to make a prediction; then (3) experimentally test the prediction. If the experiment confirms 
the prediction then declare success. If not, then go back to the drawing board. Each of these three steps 
typically requires considcrable ingenuity, insight, and resourcefulness to carry out. However, of these three 
steps, the second has a somewhat algorithmic flavour: if one only threw enough computational resources at the 
model then surely it could be solved? This echoes Hilbert's Entscheidungsproblem: perhaps for any physical 
model, when sufhciently clearly formulated, there exists a procedural algorithm to extract predictions, at least 
up to some acceptable error margin? The analogy with the Entscheidungsproblem here has been specifically 
chosen to foreshadow the answer: no, there is no such general efhcient method. 

Hamiltonian complexity has been developed to answer precisely the question: how hard is it to simulate a 
physical system ? To attack this vaguely phrased question it is convenient to introduce a decision problem to 
capture its essential elements. This problem, called the simulation problem, is described as follows. Firstly, 
we need to choose a physical system. In all cases we assume the system (usually quantum) is modelled 
with a hamiltonian H. A crucial ingredient in thc simulation problem is that we don't just have one system 
but rather the specification of a family of physical systems (usually indexed by the number, n, of 

particles involved). This is fundamental because in order to answer the question "how hard is it to simulate 
if„?" using the tools of computational complexity we are going to quantify the hardness using the scaling 
f {n) with n of the number of arithmetic operations required to simulate the system. 

Next, "simulate" is far too vague: we also need to specify what state the system is in, i. e., thermal equi- 
librium or dynamically evolving? Additionally we need to say what we wish to simulate because for systems 
comprising n particles there are typically exponentially many possible quantities which are, in principle, ob- 
scrvable, but most of these are experimentally inaccessiblc. Thus the final ingredient of the simulation problem 
is the specification of a family {A„}^2 of obscrvablcs to simulate. Finally, it is unlikely that, constrained as 
we are by thc finiteness of machine precision, we can make pcrfcct predictions for A„. Therefore we specify 
two expcctation valucs «i and a2 which our predictions must decide betwecn. With these qualifications in 
mind we now write out the simulation problem. 

Problem 3.1 (The simulation problem). 

Input: 

• A specification of a family {-ffnji^Li of hamiltonians, a specification of an initial state po;n: a- family 
of obscrvablcs {An}'^^i, a (possibly complcx) time t, two proposed expcctation valucs ai < «2, and 
an intcgcr n. 

Output: 

• "yes" if 

(3.1) (A„) < ai. 
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or 

• "no" if 

(3.2) (A„) > a2, 
where 

(3.3) (A„) = tr(p„A„), 
and 

(e^*^")Vo;„e*" 



^^•^^ ^" tr((e»*^")Vo;„e»*^")- 

Notc that. because wc allow t to bc possibly imaginary, e. g., t = i/3, the trace factor in thc denominator may 
be nontrivial. From now on wc drop the subscript n on H and p and leavc the dcpendcncc imphcit. 

Two special cases of the simulation problem havc rcccived considerable attention in the literature, namely, 
thc real-time casc whcrc t is real, and thc ground-state case where po = I/ tr(I) and p = lim^_>oo e~^^^/ tr(e~^'^^) 
which is discusscd in the ncxt scction. 

As we'll see, much cffort has bcen expcnded on understanding the simulation problem in the case where 
the observable A is the hamiltonian H of the system itself. I personally think that there is much to be 
gained by studying more general observables A: not much known about the complexity of the simulation 
problem for situations of immediate physical relevance, e. g., when A is hulk or extensive quantity such as 
thc magnctisation of a quantum magnet. But surely this should be an easy problem to simulatc? After all, 
thermodynamics tells us how bulk properties should behavc. Obviously it all depcnds on what tolerances 
e = \ai — a2\ you choose: there should be a smooth transition in the complexity of the simulation problem 
from hard to easy as e is varicd bctween, say, e = 0(2~") to e = 0{n). Indeed there is some evidence of 
such a transition: for e small cnough the problem of estimating the magnctisation is essentially equivalent in 
hardness to the ground-state problem |37| and for e ^ n, in thc case of certain regular latticc systems, there 
is an efficient approximation schcme to estimate the ground eigenvalue [14] . 

Another direction which has received considerably less attention is that of t complex, but not pure imag- 
inary, nor pure real. This problem is of considerable physical importance as it corresponds to the situation 
where a system is dynamically evolving at some finite temperature, which is the generic situation in cxper- 
iments where it is usually impossible to cool a system to its ground state. Looking further afield, it is also 
interesting to propose the study of the computational complexity of systems in the presence of decoherence. 
In both the finite-temperature and the decoherent cases there ought to be a transition in the hardness of the 
simulation problem from hard at low temperature/decoherence rate to easy for high temperature/decoherence 
rate. 



4. What is hard 

In recent years we've seen a string of countcrintuitive results presenting realistic physical systems whose 
observable properties are hard to estimate. As previously mentioned, a foundational result was presented in 
a paper of Barahona |15| who proved that for the hamiltonian H of the classical Ising model in a nonuniform 
magnetic field the ground-state (i. c, t = i lim^_>oo /3) simulation problem is NP-hard for the observable A = H 
with ai = O and «2 = 2. Barahona's construction - pertaining to a classical model which, in our quantum 
context, simply means a model whose hamiltonian is diagonal in a product basis - already presents a family 
of realistic systems whose physical properties are hard to predict: any efficient general purpose algorithm to 
simulatc the ground-state energy of Barahona's Ising model would already be capable of deciding P = NP, 
and hence doesn't exist. It is worth commcnting hcrc that these results pcrtain to the worst case where the 
model is engincered adversarily. It is fairly likely that the average case complcxity, whcrc thc interactions are 
not chosen adversarily, of a physical model will, in general, bc quitc different from thc worst casc. 

The most important foundational result for hamiltonian complexity in the quantum setting is that of Kitaev 
[57] (for a review see [7]) who showed that the ground-state simulation problem where A = H is a, 5-local 
hamiltonian (i. e., a hamiltonian for a quantum spin system whose interactions involve no more than 5 spins at 
a time) with «1 = and «2 = 0{n^'^) is complete for QMA. Kitaev's proof exploited a clever trick inspired by 
Feynman generalising an idea exploited in the proof of the Cook-Levin theorem. This construction has been 
influential in all subsequent works and can be regarded as a truly quantum generalisation of the Cook-Levin 
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theorem. Thus QMA has thc samc property as N P: any efficient classical algorithm for the ground-state 
simulation problem for Kitacv's H would yield an efRcient classical algoritlini for any other problem in QMA. 
Sincc QMA contains BQP, this would imply that all quantum computations arc classically simulablc. 

It may be reasonably argued that a 5-local hamiltonian isn't particularly natural, and this prompted 
researchers to explore exactly what kinds systems might have ground-state properties which are hard to 
approximate. Building on Kitaev's construction a sequence of results refining the original have emerged 
[S51 [Ml ESI ElEZl 13 culminating in thc surprising and countcrintuitivc rcsult |Q that therc cxist even 
one- dimensional quantum spin systems whosc ground-state energy is hard to simulatc within a prccision 
e = 0(n~"), a > 3. This rcmains truc even whcn thc system is translation invariant |36] . It is worth pausing 
a moment to appreciate the full implications of this rcsult. In thc classical casc thc one-dimcnsional ground- 
state simulation problem is in P bccause one can simply cvaluatc thc partition function via a straightforward 
product of transfer operators. The results [6] of Gottesman et. al. conclusively establish that no quantum 
analoguc of thc Standard transfer operator method cxists in thc quantum setting. 

In thc proof dcscribed in |36j it seems necessary that the local spin dimension nccds to be rather largc (at 
least more than, say, 8 or 9). It is an interesting question to study whether reducing the local spin dimension 
makes the problem easier. Some evidence this might be the case arises in the case where constraints are 
chosen randomly (see the discussion below conccrning quantum satisfiability). Here it is known |67] that 
there is a transition in the entanglement structure of the ground state of one-dimcnsional quantum systems 
as the local spin dimension changes from 2 to 4. Thc case of 3-dimensional local spins seems special. Based 
on this I would be willing to conjecture that thc ground-state simulation problem for one-dimensional chains 
of spin-1/2 degrees of freedom is in P. (Some of my coUcagucs disagrec with me that this could be truc. Any 
evidence cither way would be quitc interesting.) 

While initial research focusscd on the setting of distinguishable quantum spins, much rcmains truc in the 
physically more rcalistic setting of intcracting particle systems comprised of bosons or fermions. Along these 
lines the simulation problem for a family of systems known as stoguastic, whosc hamiltonian's ofF-diagonal 
entries are negative, has been investigated [HJ [5H [23] . Such systems do not exhibit a "sign problem" 
because the coefScients of the ground state are nonnegative. (All bosonic systems in second quantisation 
with the Standard kinetic-t-potential energy structure exhibit this property.) Similar results hold here: the 
ground-state simulation problem for stoquastic hamiltonians has been shown to be complete for the class 
stoquastic-MA which is, roughly speaking, a class generalising N P with a complexity intermediate between 
two probabilistic classes MA and AM cxtcnding NP. Additionally, thc original QMA results developed above 
have been uscd to show that ground-state simulation problem is QMA hard for both fermions |91| and bosons 

[Si- 

A very closely rclatcd problem (indecd, dual to) to the simulation problem known as thc density oper- 
ator consistency problem has also rcceivcd attention. Here the task is to decide if, given a list of reduccd 
density operators, whether there exists a global pure state consistent with these reductions. This problem 
is now known |62j to be hard for QMA. The fermionic equivalent, the N -representability problem, of central 
importance in quantum chemistry, has also been shown to be QMA hard [64]. 

Parallel to this there have been several investigations exploring physically motivated specialisations, gener- 
alisations, and alternative settings of the original local-hamiltonian results. Notable results along these lines 
include studies of thc complexity of state preparation via cooling [50] of local hamiltonian systems, and the 
difficulty of thc task of simulating the density of states [HH [IH] of a local hamiltonian. 

Therc is one special casc of the simulation problem for local hamiltonians which is of particular interest, 
namely, guantum k -satisfiability: here thc task is to decide whether the hamiltonian is frustration-free or 
not, i. e., whether the ground state is exactly thc ground state of cach interaction term, or is far from such 
a state. This situation corresponds to the simulation problem where A = H > O with «1=0 and «2 = 
l/poly(n), i. e., the task is now to simply decide if a given state is exactly a ground state or far from one 
when given the additional promise that the hamiltonian is comprised of a sum of projection operators. This 
problem is arguably a more natural quantum generalisation of 3SAT. The same proof applies, as above, for 
local hamiltonians to show that this problem is hard in exactly the same settings, i. e., for one-dimensional 
translation-invariant quantum spin systems. The case of frustration-free stoguastic hamiltonians |24j is also 
rather interesting: here one can simulate the adiabatic evolution of such models classically. 
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The quantuni A:-SAT problem has turned out to be a fruitful basis for quantum generahsations of the 
random satisfiability problem, where constraints are chosen randomly. It is now understood that so-called 
random quantum /c-SAT exhibits a transition [SU [60l EH EH El] , analogous to its classical counterpart [56] , 
from a phase where there are almost always many possible frustration-free ground states to a phase where 
there are almost no frustration-free ground states. 

A setting which has been the subject of considerable conjecture and interest is that of quantum spin 
Systems with commuting interactions |251 [S] [B^. One might be forgiven for initially guessing that such 
Systems are completely classical because all interactions can be simultaneously diagonalised. However, a little 
thought shows that this simultaneous diagonalisation isn't particularly hclpful because the basis in which all 
the interactions are trivial may be highly entangled. Thus, the ground-statc simulation problem may well be 
hardcr than its classical counterpart. Indecd, some evidence this might be the case has been recently obtained 
[5]. Deciding whethcr the simulation problem for commuting interactions is harder than NP, or in NP, is an 
intriguing opcn problem because it is clear that it rcquires new ideas going beyond the existing tcchniques. 

4.1. Guiding problem: the quantum PCP conjecture. To date, in all of the constructions of quantum 
Systems whose ground-state energy is hard to simulate, the precision e = |qi — 02], or promise gap, for which 
the problem is known to be hard scales as an inverse polynomial in the system's size. That is, the problem is 
known to be hard only when one requires that the energy is very precisely approximated. It is thus natural 
to ask to what extent the intractability of these systems is tied to this constraint, i. e., perhaps it is easy to 
decide if the ground-state energy is small or larger than, say, some constant fraction of the number m of 
interaction terms? 

In the classical setting this question has a dramatic and counterintuitive answer: it is also complete for N P. 
This result, captured by the celebratcd PCP theorenjl [EKTTIIST], shows that even the task of approximating 
the ground-state energy to within a constant factor is hard. 

The PCP theorem already establishes that there are physical systems, indeed classical systems, for whom 
the simulation problem is hard, even when e = \ai — «2! = cm, with c some constant. However, it is entirely 
natural to suppose that that there should be a corresponding quantum equivalent of the PCP theorem which 
would say that the ground-state simulation problem is QMA-complete for e = jai — 02] = cm, with c some 
constant [31 [5]. 

Some progress towards a possible proof of the conjecturcd guantum PCP theorem has been made whereby 
an important proof element known as gap amplification - exploited in a recent proof [31] - has been generalised 
to the quantum setting, where it is now known as the detectability lemma [3]. However, so far, a fully quantum 
PCP theorem has resisted all attempts at a proof. 

The quantum PCP conjecture is a fascinating guiding problem because it would require the development 
of many theoretical tools of indcpendent interest. For example, it is considered likely by practitioners that 
a proof of the quantum PCP conjecture would aid in the design of error-tolerant architectures for adiabatic 
quantum computation. Additionally, proof techniques developed in the pursuit of the quantum PCP con- 
jecture also have direct applications in the study of the local hamiltonians appearing in condensed matter 
physics. For example the detectability lemma has already found applications in the study of entropy/area 
laws for quantum spin systems (see below) |3]. A quantum PCP theorem would mean that there exist certain 
realistic hamiltonians for which no variational representation of the ground state with a concise classical 
description could reproduce the ground-state energy bettcr than a constant precision. 

At this point it seems that a proof of the quantum PCP conjecture would have an inductive structure 
whereby a very bad quantum PCP is iteratively improved. However, there is a stumbling block in any 
proof following this strategy, namely, the design of at least one, possibly very inefficient, quantum PCP. 
This is summarised by the following problem: show that a QMA-complete simulation problem for n spins 
can be encoded as the simulation problem for a local hamiltonian involving an exponential (or even doubly 
exponential) number N of spins, i. e., N = 2", oi N = 2^ , with e = \ai — «2] — cN, c a constant. In other 
words, if we don't carc how many spins we use is it possible to improve the promise gap e ~ \ai — «2! at all 
while retaining the hardness of the problem? Historically this was one of the stcps toward the classical PCP 
theorem. What I have in mind here is something like a quantum generalisation of the results of Hastad [10] 
on inapproximability of certain problems involving, e. g., solutions to linear equations. 



•^PCP stands for probabilistically checkable proof. 
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5. What IS EASY 

In the previous section we saw tliat considerable efFort has been expended on the task of constructing 
plausible physical models whose ground-state properties are hard to approximatc. In this section I'll takc a 
look at the flip side of the coin, i. e., understanding what physical properties are easy to approximate. To 
answer this question concretely requires the description of an algorithm which is correct, converges, and has 
a polynomial complexity. The design of such algorithms is truly an art form, but one informed by physical 
results. 

The simple and striking conclusion we've drawn so far is that the natural dynamics of a locally interacting 
System doesn't always take the system an equilibriuni state modelled by the Gibbs distribution. The question 
then arises: what states can a locally interacting system enter under natural dynamics? The key insight to 
draw on here is that the dynamics of a locally interacting coUection of particles cannot explore the entirety 
of hilbert space in any reasonable length of time as it takes only a polynomial number, in n, of parameters 
to completely specify a local hamiltonian, but to specify a state in hilbert space requires an exponential 
number of parameters. Thereforc, in order for a local system to cxplore hilbert space via its local dynamics, 
it ought to require at least an exponcntial time to visit every location. (Here we are thinking of hilbert space 
as a manifold rather than a linear space, i. e., we are forgetting the linear structure.) This simple counting 
argument can be improved |83| and it turns out that the strongest result is much better: it takes a doubly 
exponential time for the dynamics of a local hamiltonian to visit an arbitrary state in hilbert space. 

To specify a pure state in the classical situation one requires only a linear (in n) number of parameters, 
namely, the positions of each of the particles. Another way of saying this is that for any pure classical state 
one can assign a dcfinite pure state to each of the constituent particles. The main difhculty one faces in 
simulating a many particle quantum system is that it is difficult to describe a generic pure state in the hilbert 
space of n particles, because a general such state can be a superposition of possibly exponentially many 
classical configurations. I. e., for a general pure quantum state one cannot assign a definite pure state to each 
of the constituent particles as the parts do not completely determinc the whole. The mismatch here is due 
to guantum entanglement, coming from the quantuni corrclations between the constituent particles. Thus, in 
order to completely describe the reduced state of a single particle we also nced to keep track of the quantum 
corrclations that devclop between it and the rcst of the system. This is, generically, difficult. 

Howevcr, we can cxploit locality structure to placc bounds on how fast Information can propagate through 
the system under natural dynamics. The physical picture to have in mind here is that of a lattice of particles 
connected by springs: if you wiggle one particle it should take a time proportional to the distance multiplied 
by a "speed of sound" for your disturbance to propagate through the system. This picture does remain 
metaphorically correct in the quantum setting and a fundamental result known as the Lieb-Robinson bound 
quantifies this intuition: suppose that H is the hamiltonian for a quantum spin system on a regular lattice 
and A and B are two observables acting nontrivially on two disjoint regions, then [ST] (see [T7] for a textbook 
treatmcnt) 

(5.1) \\[e''"Ae-''",B]\\ < c|| A|| ||S||e"l*l-«^, 

where c, w, and k are constants, and L is the distanc^^ between the regions that A and B act on. The Lieb- 
Robinson bound says that regardless of the initial state of the system, if Alice implements some operation at 
time r = O then no matter what observable Bob subsequently measures at some later time T = t the amount 
of Information that Alice can send to Bob is exponentially suppressed as the separation between them is 
increased [IH]. Note, however, that as t increases the bound becomes cxponentially weakcr. This trade-off 
between exponentials gives rise to an effective "Information propagation cone" : Information about Alice's 
operation is exponentially suppressed outside of the region {x\d{x^A) < -\t\}, where d{x,A) measures the 
distance between site x and Alice's region. 

The proof of Lieb-Robinson bound has been simplified and generalised to a variety of scttings beyond the 
original translation invariant setting (see, e. g., |73| for a recent review): beginning with |65[l86j we've witnessed 
a sequence of results taking [HJ H?! [TU |BS1 [53] it beyond the translation-invariant setting to systems with 
disorder [27] , classical systems [87l [9] , unbounded operators [28l [30] [69] [TO] [84] [85] , and dissipative systems 



Measured using the natural lattice distance. 
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[82l [74] . Applications of thc Lieb-Robinson idea have had influences further afield, e. g., to givc ncw proof 
techniques for understanding the emergence of effective equations in mean-field settings [34] . 

Armcd with the Licb-Robinson bound we ncw havc the necessary quantitative tool to undcrstand cxactly 
how to parametrise naturaUy occurring states. Intuitively, although each particle does beconie correlated 
with the others after a smah time, it should only be able to develop correlations with its neighbours within 
a distance proportional to the Lieb-Robinson velocity. One can then keep track of only these correlations 
and obtain a compact description of the global state. Thus we obtain an ejfectively local dcscription of the 
system's state, i. e., we speeify the state by reeording the correlations between any given partiele and only its 
immediate neighbours. The remarkable eonsequence |78[|44[H5] is that, in one dimension, there is an efficient 
representation known as a matrix product state (MPS) for the system's state, as long as |t| < log{N). An 
algorithm exploiting this phenomenon was already proposed in [97] : the preceding results ean be interpreted 
as saying that this algorithm (and relatives) is indeed correct, converges, and has polynomial complexity. 
That is, the real-time simulation problem is in P. Further advances along these lines have shown that the 
simulation problem for adiabatic evolution of gapped quantum spin systems also lies within P [79]. 

The Lieb-Robinson teehnology has also been put to great effeet by Hastings and eoworkers in the study of 
the equilibrium properties of quantum lattice systems with a spectral gap between the lowest two eigenvalues. 
Beginning with the breakthrough result [41] , where the the Lieb-Schultz-Mattis conjeeture was proven in higher 
dimensions, mueh insight has been gained into the general properties of gapped quantum lattiee systems (see 
also [75]). Firstly, the presence of the gap allows one to play the Lieb-Robinson bound off against the natural 
decay induced by the gap to show that all spatial ground-state correlators deeay exponentially with separation 
[4T1 [42l [TTI [48[ [39[ [13] . This result already hints at a more powerful generalisation: perhaps the entanglement 
entropy, as measured by the von Neumann entropy S{pl) ~ — tr(pi log(pL)), where is the reduced density 
operator of a region of length L in the ground state of a gapped spin system, is bounded by a constant? 
Indeed this is the case for one-dimensional systems and sueh an entropy /area law was proved in [43]. This 
result can be interpreted as placing the ground-state simulation problem for gapped ID spin systems in N P, 
i. e., there exists a polynomial sized proof, written as an MPS, that the ground-state energy is low. This 
problem is not known to be in P. 

The proof technique developed in [43] also established the key result that the ground-state of a gapped one- 
dimensional quantum lattiee system may be efficiently approximated by an MPS. The proof is not eonstructive, 
however, and only supplies a promise that such a representation exists. It is tempting to suppose that a 
simple application of the variational method to MPS should find such a representation easily. However, a 
naive application of the variational prineiple to MPS runs into problems bccause the objective function is 
a nonlinear funetion with many loeal minima, and it is far from clear how to carry out the minimisation 
efficiently. 

It turns out. however, that a practical solution to this problem has been known for many years under 
the guise of the density matrix renormalisation group (DMRG) [99| . The DMRG is now understood to be a 
relazation of the variational method applied to MPS [88]. While the method is a heuristie, and it does sufFer 
from loeal minima problems, extensive praetieal numerical experience has shown that these local minima 
do not cause any significant problems. Indeed, the DMRG is, without peer, the premiere method for the 
exploration of strongly correlated quantum physics in one dimension and is hard to overstate the impact of 
the DMRG in thc field of condensed matter physies. 

The question of whether there does exist a polynomial algorithm to vary over MPS has, however, been 
addressed. It was recently established [2] [90] that a transfer-operator type argument yields an efficient 
method to vary over MPS with a fixed number of variational parameters. Unfortunately, the performance of 
this algorithm is not eomparable with the DMRG, and is unlikely to supplant the DMRG as thc method of 
choice in the study of one-dimensional quantum systems. 

The MPS variational class was proposed and studied well before the advent of hamiltonian complexity. 
But with the advent of quantum Information theory and hamiltonian eomplexity new ideas began to flow into 
condensed matter physies leading to the development of several new variational elasses for quantum lattiee 
systems. One important example are the projected entangled-pair states (PEPS) |93j . whieh are a natural 
generalisation of MPS to the two-dimensional setting. Another example is the multiscale entanglement renor- 
malisation ansatz (MERA) [M] [M] , whieh is particularly well-suited to the description of low-dimensional 
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quantuni lattice systems at or near cricality. Both of these classes have been put to great cffcct, in combination 
with thc variational principle, in the practical study of quantum lattice systems. 

Parallcl to thcsc theoretical developments thcrc has been the exploration of classical heuristics and algo- 
rithms to approxiniate ground-state properties of quantum spin systems. Firstly, it was reahsed |20| that one 
can exploit the low local dimensionality of quantum spins to efficiently decide the quantum 2-satisfiabihty 
problem in the spin-1/2 case. A variation on this theme exploits possible low-degree structure an interaction 
hypergraph to supply a criteria to decide whether a quantum satisfiability instance always has a solution. This 
result [5] is known as the guantum Lovasz local lemma. Further afield, the study of approximation algorithms 
and heuristics for QMA hard instances of the local hamiltonian problem has been inccpted in I35| . 

Recently, in another direction, Lieb-Robinson ideas and a technique known as quasi-adiabatic continuation 
[47l [T3] have been exploited to prove several powerful results concerning the stability of topological order 
[22l [18] . Quasi-adiabatic continuation has been further explored in the contcxt of disordered systems [46] , 
where the improved bound on Information flow coming from the Anderson localisation phenomena may 
be exploited to improve many results concerning the decay of correlations in strongly interacting quantum 
systems. One of the most recent triumphs of these research directions has been the proof [32] of an outstanding 
conjecture in mathematical physics, namely, the quantisation of the Hall conductance without averaging 
assumptions. 

5.1. Guiding problem: entanglement entropy/area laws. The entanglement entropy/area law states 
that for a quantum system in the ground state the entropy of any subregion A should scale as the area of the 
boundary. This result may be regarded as an expression of the holographic principle, namely, that the physics 
in the bulk of the system can be completely described by what happens at the boundary. The existing proofs of 
these results have been influential in hamiltonian complexity and condensed matter physics, typically leading 
to many new insights and results, including new efhcient methods to simulate strongly correlated systems. 

Quantitative results concerning the entropy/area law are, however, difficult to obtain. Indeed, apart from 
special cases |66j . the only settings where the result has been proved at the level of mathematical rigour 
pertain to gapped one-dimensional quantum spin systems and gaussian systems (see, e. g., [32] for a review). 
One generalisation that has proved particularly elusive is to gapped two-dimensional quantum spin systems. 
Here the existing proof ideas are extremely suggestive, yet just fail to yield the expected law. The difhculty is 
probably tied, in part, to the fact that two-dimensional quantum systems can exhibit topological order, which 
is difficult to capture using the existing proof techniques. It is clear that whatever method is developed to 
tackle this case will have important consequences not only for the simulation of two-dimensional quantum 
systems, but also for the study of exotic phases of matter. 

6. CONCLUSION 

Hamiltonian complexity, as a field, has developed rapidly in the past couple of years amassing a broad 
coUection of foundational results exploiting a wide portfolio of techniques. Highlights include the proof 
that it is QMA-complete to approximate the ground-state energy of even a one-dimensional quantum lattice 
system and, in the other direction, proofs that the imaginary time simulation problem is in N P for gapped 
one-dimensional quantum lattice systems, and the real-time simulation problem is in P for \t\ ^ log(n). 

The task of finding quantum lattice systems which are hard to simulate appears to have reached a natural 
pause, of sorts. This is probably due to the difRculty of the proofs required and the fact that the existing 
techniques have been pretty much pushed as far as they can go. However, several simple-sounding problems 
remain, such as the commuting hamiltonian problem, and the search for inefficient quantum PCPs, and 
progress on these would open up many new directions. 

The development of efficient algorithms for thc simulation of low-dimensional quantum lattice systems 
continucs unabated. In this direction I think the study of quantum approximation algorithms initiated in |35| 
holds particular promise. 

One area that I am particularly excited about is that of the generalisation of thc insights obtained in 
hamiltonian complexity to quantum systems with continuous degrees of freedom, including, quantum fields; 
we've recently understood [M] (HOl [3H] how to generalise the two main variational classes studied in the 
hamiltonian complexity literature, namely MPS and MERA, to the continuous setting, and it seems likely 
that many new results generalising those originating in the lattice setting will be possible. 
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